Dimension-Dependence of the Critical Exponent in Spherically Symmetric 

Gravitational Collapse 



in 
o 
o 

(N 

< 
O 

m 

> 
00 
00 

o 
l> 
o 
in 
o 

cr 

i ' 



X 



J. Bland" and B. Preston*, M. Becker", G. Kunstatter", V. Husain b 
bDept. of Mathematics and Statistics 
University of New Brunswick, 
Fredericton, N.B. Canada E3B 1S5. and 
ftDept. of Physics and Winnipeg Institute of Theoretical Physics 
University of Winnipeg, Winnipeg, Manitoba Canada R3B 2E9. 
(Dated: February 7, 2008) 

We study the critical behaviour of spherically symmetric scalar field collapse to black holes in 
spacetime dimensions other than four. We obtain reliable values for the scaling exponent in the 
supercritical region for dimensions in the range 3.5 < D < 14. The critical exponent increases 
monotonically to an asymptotic value at large D of 7 ~ 0.466. The data is well fit by a simple 
exponential of the form: 7 ~ 0.466(1 - e -0 ' 40813 ). 

PACS numbers: 04.70.Dy 

I. INTRODUCTION 

ChoptuikQ was first to discover that gravitational col- 
lapse to black holes exhibits critical behaviour 0, 0]. In 
particular, if one evolves initial data corresponding to a 
matter distribution described by some parameter a, there 
is a critical value a» of that parameter that separates data 
whose final state is a black hole, from that for which the 
matter ultimately disperses to infinity. Near the critical 
value, the radius of the horizon on formation obeys the 
simple scaling lawpj: 



Rah(a) = (a - a*)"- 



(1) 



part by the relationship between critical phenomena in 
gravitational collapse and black hole to black string phase 
transitions in higher dimensions, Sorkin and Oren[9( per- 
formed a similar analysis using a different formalism. 

In the following, we describe the dependence of the 
critical exponent 7 on spacetime dimension. A simple, 
but significant, improvement to the original code in [8j 
has enabled us to obtain reliable and accurate results for 
dimensions 3.5 < D < 14 [l^. In this range the critical 
exponent grows monotonically with D to an asymptotic 
value near 0.466. Remarkably, the data is well fit by a 
simple exponential: 



The scaling exponent 7 depends only on the type of mat- 
ter that is considered (eg scalar field, vector field, etc.) 
and not on the choice of parameter a or the distribu- 
tion of the initial data. The emergence of this universal 
exponent is fairly well understood numerically in terms 
of a unique, discretely self-similar critical solution which 
serves as an attractor. The self-similarity is manifested 
by a periodicity of the critical solution for the scalar field 
x(r, t) of the form: 



7 = 0.466(1 - e 



-0.4081) 



(3) 



X(r,t) = x(re A ,te A ) 



(2) 



with echoing period A. The critical solution has only 
been studied numerically and very little is known analyt- 
ically about it, except perhaps in the context of the three 
dimensional BTZ black holejj]. 

Until recently, little was known about critical grav- 
itational collapse in spacetime dimensions other than 
four. The only exceptions were in three^, 0] and six 
dimensions^- A few years ago, two of the current au- 
thors started a program [8( designed to examine the de- 
pendence of the critical exponent on spacetime dimen- 
sion. A formalism was set up that allowed the use of 
a single code in which dimension appeared as an input 
parameter. Previous results in four and six dimensions 
were confirmed and some preliminary new results in five 
dimensions were presented. More recently, motivated in 



Our data are consistent with those of 0. Our error 
bars are however smaller, and we are able to get reliable 
results in higher dimensions. We see no evidence for the 
maximum in 7 near D — 11 suggested in [ij. 

We first review our general formalism, then describe 
the numerical method and calculation of error. Finally, 
we present our results and conclusions. 



II. FORMALISM 

We begin with Einstein gravity in D spacetime dimen- 
sions minimally coupled to a massless scalar field: 
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(4) 



Spherical symmetry is imposed by writing the metric 
g/ii' as 

ds2 (D) = g a pdx a dx 13 + r 2 (x a )dn {D - 2 ), (5) 
where dClm-2) is the metric on S D ~ 2 and a, j3 = 1, 2. 
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The key to obtaining a simple and elegant set of equa- 
tions in all dimensions is a field redefinition motivated by 
2-d dilaton gravity [j [TQj: 



9aP 



8(n- 1) 

.(n-l)/n 



9af3, 



(6) 
(7) 



where n = D — 2. Note that r is the optical scalar, and 
is proportional to the area of the n-sphere at fixed radius 
r. Formation of an apparent horizon in such a spherically 
symmetric spacetime is signalled by the vanishing of the 
quantity: 



ih = g a ^d a 



(8) 



We go to double null coordinates first introduced by 
Garfinklejl^: 

ds 2 = —2lg(u,v)(f>'(u,v)dudv (9) 
in which the relevant field equations take the simple form: 

I 

= o*x\X ) 

(12) 



0' = _i.yW(^' 
WY + (4>x)' = o, 



(10) 

(ii) 



where prime and dot denote the v and u derivatives re- 
spectively. All the information about spacetime dimen- 
sion is encoded in the potential VW; 



with: 



vWfaA) = (n-1) 



8(n- 1) 



l 2 Acj) 



l/r, 



(l)' 



(13) 



(14) 



As one might expect from the fact that critical collapse 
is a short distance phenomena, the cosmological constant 
does not affect the critical exponent. This was verified 
by explicit calculation in We will henceforth set 

A to zero and concentrate on the dependence of 7 on 
spacetime dimension D. 

The evolution equations may be put in a more useful 
form by defining the variable 



20 X ' 



(15) 



which replaces the scalar field \ by h. The evolution 
equations become 



-9/2 



(16) 
(17) 



g = exp 



47T / dvUh~X? 



X 



i\R> Ju 



dv 



h<p' 



(18) 
(19) 

(20) 



This was the final form of the equations used for nu- 
merical evolution in 8]. Note, however, that the inte- 
grand of g is singular in the critical region near tf> = 0, 
even though the integral itself is well behaved. This leads 
to instabilities that severely limit the accuracy of the cal- 
culation. A drastic improvement can be obtained by ob- 
serving that: 



V gV^U'dv = (n-1) ( -^-) 



— l/n r v 



-1/n, 



8n 
n-1 



— l/n „ v 



1-1/n 



(21) 



so that a straightforward integration by parts yields: 



9 = 



1<j)V { 



n) 



Aim 



gV^(h- X ) 2 4>'dv (22) 



n — 1 n — 1 

where we have used the expression for g in (|18f) to write: 



dg = Ang(h - X f 



-dv 



(23) 



The first term in ill'L'l) vanishes at <f> = 0, and the inte- 
grand in the second term now vanishes at the origin by 
virtue of JT^J (assuming <j>' ^ 0). Thus the use of 
instead of (|19f) improves the stability of the code near 
criticality significantly and allows reliable results to be 
obtained in dimensions up to at least 14. 



III. NUMERICAL METHOD 

The numerical scheme uses a v ('space') discretization 
to obtain a set of coupled ODEs: 



h(u, v) — > hi(u), (j)(u, v) — > 4>i(u). 



(24) 



where i — 0, • • • , N specifies the v grid. Initial data for 
these two functions is prescribed on a constant v slice, 
from which the functions g(u, v), g{u 1 v) are constructed. 
Evolution in the 'time' variable u is performed using the 
4 th order Runge-Kutta method. The general scheme is 
similar to that used in |l4j, together with some refine- 
ments used in [l^. This procedure was also used for the 
3— dimensional collapse calculations in 6J. The grid size 
used for all calculations reported here was 6000, with no 
mesh refinement. An improvement in accuracy was, how- 
ever, obtained by adaptively decreasing the time steps 
once apparent horizon formation commenced, i.e. once 
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the function ah(<p) developed an extremum. Time steps 
ranged from Au — 10~ 3 at the beginning of the calcula- 
tion to about Au = 10~ 9 near the end. 
The boundary conditions at fixed u are 



D = 4 (outer minimum) 



* = 0, h = 0, g k = 1. 



(25) 



where k is the index corresponding to the position of 
the origin = 0. All grid points < i < k — 1 corre- 
spond to ingoing rays that have reached the origin and 
are dropped from the grid. The boundary conditions are 
equivalent to r(u,u) — 0, g\ r =o — <?(m, u) = 1. Notice 
that for our initial data, 4>k and hence hk, are initially 
zero, and therefore remain zero at the origin because of 
p7jl. The initial scalar field configuration x{ u — 0, <j>) is 
most conveniently specified as a function of <fi rather than 
r. (Recall that <f> oc r n .) This together with the initial 
arrangement of the radial points 0(0, v) fixes all other 
functions. We used the initial specification <fi(0,v) = v. 
We considered Gaussian initial data 



Xg(u = 0,0) = a0 



2/n 



exp 



(26) 



The initial values of the other functions were determined 
in terms of the above by computing the integrals for g n 
and g n using Euler's method for equally spaced points. 

At each u step, a check is made to see if an appar- 
ent horizon has formed by observing the function ah in 
© whose vanishing signals the formation of an apparent 
horizon. For each run of the code with fixed parameters 
a and a, this function is scanned from larger to smaller 
radial values after each Runge-Kutta iteration for the 
presence of an apparent horizon. In the subcritical case, 
it is expected that all the radial grid points reach zero 
without detection of an apparent horizon. This is the 
signal of pulse reflection. In the supercritical case, black 
hole formation is signalled by ah vanishing at some finite 
radius. Since the code crashes before ah — is actually 
reached, one has to terminate it when ah(tfi m i n ) is below 
some bound. Note that by virtue of our parameteriza- 
tion, in vacuum ah decreases monotonically to zero at 
r = 0. As the radius of horizon formation decreases, it is 
therefore necessary to decrease the value of the bound. 
The required values ranged from about 10 -6 near criti- 
cality in lower dimensions to 10~ 34 in higher dimensions. 

In all cases, we used values of 0o = 1 and a = 0.3 for 
the Gaussian initial data and varied the amplitude a to 
determine the critical amplitude o* below which no black 
hole formation was observed. We then did runs in the 
supercritical region a > a* to determine the dependence 
of the horizon radius R a h on the amplitude. The critical 
exponent is given by the slope of the best fit straight line 
for the hi(R a h) vs ln(a— a*) plot, while the echoing period 
is determined from the periodicity of the residuals. See 
Fig.®. 

Universality was tested by fixing a = a* in the pro- 
gram, treating a — er» = 0.3 as the critical width, and 
then varying the width in the supercritical region a < a* . 
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FIG. 1: Plot of In(RAH) vs. ln(a — a,) for 4 dimensional 
case. The residuals clearly show the oscillation expected due 
to discrete self-similarity of the critical solution. 



In all cases that we tested, varying the width and vary- 
ing the amplitude gave the same critical exponents to the 
required accuracy. For example in 12 dimensions (see 
Fig-©) the difference between the two values was less 
than 0.5%. 



D = 1 2 (width) 




FIG. 2: Plot of Iti(Rah) vs. ln(a — a,) and ln(R a h) vs. 
ln(a, — a) for our 12 dimensional case. Although the residuals 
are cleaner in the amplitude case, the slopes of the lines agree 
to within error. 



The accuracy of the critical exponent is in principle 
limited by the number of significant digits in the critical 
amplitude. Due to instabilities in the code near critical- 
ity, it became progressively more difficult to determine 
a* directly as the dimension increased. Our error esti- 
mates were determined as follows. We first determined 
a range for a* by direct search. We then varied a* in 
the determined range, and for each value we calculated 
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and plotted Iu(Rah) vs ln(a — a*). The value of a* used 
to determine 7 as quoted in Table (1) was the one for 
which the average squared deviation R 2 of the Iu(Rah) 
vs ln(a — a») plot was minimized. Since the tail of the 
Itl{Rah) vs ln(a — a*) plot is fairly sensitive to the value 
of a* , we were able to narrow the range of a* , and hence 
lower the error bar on the slope, by determining the val- 
ues of a* for which R 2 was within a specified range of 
its maximum (generally ±0.0001). Although there was 
some variation with dimension, the procedure generally 
gave us an error in 7 of about 1%. The error was less in 
4 dimensions, for which the code was most stable. It is 
interesting that the overall error did not grow rapidly at 
high dimension, despite the fact that our determination 
of a* became less accurate. 

It is also important to note, that in all cases tested, 
including 14 dimensions, varying the width of the initial 
matter distribution yielded the same value for the critical 
exponent as varying the amplitude, within the given error 
bars. See for example Fig.J2J. 




FIG. 4: Plot of ah vs. phi for 10 dimensions, showing the ap- 
pearance of multiple minima. The outermost minimum cor- 
responds to the horizon location. 
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FIG. 3: Plot of Iu(Rah) vs. ln(a — a*) for D = 10. The 
curve on the left tracks the radius of the true horizon on 
formation, whereas the curve on the right tracks the same 
minimum throughout. The slopes are within error, but the 
residuals on the right clearly give the echoing period associ- 
ated with discrete self-similarity. 



Although the horizon location is always given by the 
outermost minimum the recorded value jumps as new 
minima appear. When the location of a single minimum 
(not necessarily the outer one which formed the horizon) 
is tracked throughout the range of amplitudes considered, 
the residuals are much smoother, and give a more accu- 
rate indication of the echoing period. This procedure 
does not significantly affect the slope (see Fig.©). 

Another, potentially more accurate, estimate for A can 
be obtained by observing the ringing of the scalar field 
near the origin just prior to the formation of an apparent 
horizon near criticalitv|l7j. Although our graphs gave 
a clear indication of the ringing (see Fig. (jSJ), we were 
unable to get close enough to criticality to obtain very 
accurate results for the period. This quantity is generally 
easier to obtain by observing the scalar curvature in the 
sub-critical region, as done in 9]. We estimate that our 
error in determining A was in the neighbourhood of ±0.1. 



IV. RESULTS AND DISCUSSION 



The echoing period can be obtained in a straightfor- 
ward way by measuring the period of oscillations of the 
residuals from the straight line in the ln-ln plots. It is 
given by: 

A = 2T 7 (27) 

where T is the oscillation period. The accuracy is lim- 
ited by the errors in both 7 and T , and is therefore 
considerably worse than that of the critical exponent. 
The determination of A is further complicated in higher 
dimensions by the appearance of multiple minima (see 
Fig.©). 



The results for the critical exponent and echoing period 
are given in Table (1) below. The universality of the 
results, plus the clear presence of oscillations in the ln-ln 
plots confirm that we are sufficiently close to criticality 
in all cases to obtain reliable results. 

The critical exponent appears to grow with dimension 
monotonically to an asymptotic value near 0.466. As 
seen in Fig. 10, the data are well fit by an exponential of 
the form: 7 = A - Be~ Cy D with A ~ B ~ 0.466 and 
C ~ 0.408. The fact that A and B are approximately 
equal is significant, since they are a priori independent 
parameters. 

Our results are consistent within error with those ob- 
tained by Sorkin and Oren[9j, for both the critical expo- 
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X Near Origin in 13 Dimensions 



Dimension Dependence of Scaling in Spherically Symmetric Scalar Field Collapse 
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u, = 2.4856 + 0.0006 



2 3 4 5 

-ln(u. - u) 




4 6 8 10 12 

Spacetime Dimension, D 



y= 0.466(1 - exp[-0.408D]) 

14 16 



FIG. 5: Plot of x vs - m (« — u *) near criticality for 13 dimen- 
sions, showing the ringing expected from the discrete self- 
similarity of the solution. 



FIG. 6: Plot of results for 7 as a function of spacetime di- 
mension. 
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TABLE I: Results for critical exponent and echoing period. 
Universality of the critical exponent was checked in several 
dimensions by varying width and amplitude. Comparisons 
with previous work of |j is given in columns 3 and 5. 



nent 7, and the self-similar period A. However, we do 
not observe the maximum for the critical exponent near 
D = 11 reported in There is no contradiction here, 
since the authors of that paper calculated the critical 
exponent to 11 dimensions, and conjectured a maximum 
by extrapolating polynomial fits to higher dimensions. In 
this context it may also be useful to note that a maximum 
in the critical exponent was also present in the prelimi- 
nary results reported by one of the authors in jl .'J| . In our 
case, the maximum was determined to be an artefact of 
the inaccuracy of the calculation. At higher dimensions 
it is more difficult to obtain numerical stability near crit- 
icality. Moreover, the slope of the plot of ln(R ah ) vs. 
ln(a — a*) systematically decreases as you move further 
away from a*. The critical exponents reported in |l5j | de- 



creased for dimensions greater than ten because the data 
were generated by points progressively further away from 
criticality. 

Our results for the ringing period A suggest that A de- 
creases monotonically to an asymptotic value of just be- 
low 2.4. This is not unreasonable given our more reliable 
results for the critical exponent. If indeed the D = 00 
limit of the questions is well defined, then it is reasonable 
that both 7 and A asymptote to a finite value determined 
by the critical solution of the limiting theory. It is of 
great interest to try to find this solution analytically. In 
fact one of the original goals that motivated the program 
started in Q was the hope that a numerical analysis of 
the dimension dependence of critical collapse would pro- 
vide clues about the nature of this limiting solution. This 
is currently under further investigation. 

In summary, we have presented the results of an anal- 
ysis of the dependence of the critical exponent in spher- 
ically symmetric scalar field collapse on spacetime di- 
mension. There is strong evidence that this exponent 
increases monotonically to an asymptote around 0.466. 
The data are remarkably well fit by the simple relation- 
ship 0, which also suggests that the critical exponent 
is finite at D = 3. However, one can also fit the data, 
albeit not quite as well, (see Fig.JTJl), with an expression 
of the form: 

7 = 0.4727- (28) 

which goes to minus infinity at D = 3. In order to set- 
tle this point definitively we need to explore the region 
near D = 3 more thoroughly. This work is currently in 
progress.. 

Given the recent interest in higher dimensions in gen- 
eral, and higher dimensional black holes in particular it is 
worthwhile to understand in detail the nature of gravita- 
tional collapse to black holes in all dimensions. This may 
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FIG. 7: Alternative fit for 7 as a function of spacetime di- 
mension that diverges at D = 3. 



be particularly relevant in the context of black string to 
black hole transitions in higher dimensions. The results 
presented here should provide valuable clues in address- 
ing these issues. 
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